Pulse propagation in tapered granular chains: An analytic study 
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We study pulse propagation in one-dimensional tapered chains of spherical granules. Analytic 
results for the pulse velocity and other pulse features are obtained using a binary collision approxi- 
mation. Comparisons with numerical results show that the binary collision approximation provides 
quantitatively accurate analytic results for these chains. 
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I. INTRODUCTION 

The study of pulse propagation through granular me- 
dia has attracted a great deal of attention for several 
reasons. Apart from addressing fundamental problems 
of pulse dynamics in the presence of highly nonlinear in- 
teractions, it also has direct practical application. Pulse 
propagation has mostly been studied in monodispersc 
one-dimensional chains of spherical granules. Nesterenko 
first showed P, that an initial impulse at an edge of a 
granular chain in the absence of prccomprcssion can re- 
sult in solitary waves propagating through the medium. 
Since then, this system and variants thereof have been a 
testbed of extensive theoretical and experimental stud- 
ies [I lip, S 0, S i, M, ini, E [H, M M: M, [13, [H, 

m, 0,i2ir. 

The chain model has recently been generalized to study 
the (sometimesprofound) effects on pulse propagation of 
dissipation and of polydis pers ity in the struc- 

ture and mass of the granules [1, QlH, [S, [13, [H, [H. 
In particular, polydispersity is frequently introduced in 
a regular fashion such as in tapered chains (TCs), in 
which the size and/or mass of successive granules sys- 
tematically decreases or increases. Polydispersity is also 
introduced by distributing masses randomly, by "deco- 
rating" chains with small masses regularly or randomly 
placed among larger masses, and by optimizing grain dis- 
tribution for particular purposes. For example, consider- 
able recent activity has focused on pulse propagation in 
"forward tapered chains" [H, [H, [ij, Il9l . with some 
results corroborated experimentally [20, l2l|. In these 
chains each granule is smaller than the preceding granule 
according to a systematic pattern. Among the main re- 
sults of these studies is that the energy imparted to the 
first granule is redistributed among increasingly larger 
numbers of granules as the pulse propagates, thus atten- 
uating the impulse on subsequent granules as time pro- 



ceeds. This is a desirable behavior in the context of shock 
absorption. In any case, almost all known results are nu- 
merical and therefore predictively limited. While a few 
specific realizations have been experimentally verified, we 
have found few analytic results for pulse propagation in 
polydisperse systems [isj . Analytic work is in principle 
possible (albeit in practice often extremely difficult) if a 
continuum approximation is valid. This approach was 
first implemented successfully for monodispersc granular 
chains, where the resulting continuum equations can be 
(approximately) solved analytically [J, [3, H, H, 0] . Even 
when a continuum formulation is appropriate, in many 
cases the resulting equations can not be solved [1, ■ 

In this work we develop a converse approach which 
focuses on the fact that the propagating front in granu- 
lar chains is often narrow. Our theory invokes a binary 
collision approximation which leads to analytic expres- 
sions for various quantities that characterize the prop- 
agating pulse. In an earlier paper we introduced a bi- 
nary collision approximation for the calculation of the 
pulse velocity in monodispersc chains Interestingly, 
although it is "opposite" to the continuum approxima- 
tion, we found that the binary collision model produces 
results that also capture the dynamics of pulse propaga- 
tion in such chains. Here we generalize this approach to 
forward and backward TCs, that is, to chains in which 
successive granules increase in size or decrease in size. 
The high quality of the approximation is ascertained by 
comparison with numerical integration of the full equa- 
tions of motion for two tapering protocols. In one, the 
granular radius increases or decreases linearly (backward 
or forward "linearly" TC), while in the second the ra- 
dius increases or decreases exponentially (backward or 
forward "exponentially" TC). 

In Sec. In] we introduce the granular chain model 
and lay out the binary approximation procedure and its 
caveats. The analytic results for the backward TCs and 
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comparisons for these chains with numerical integration 
results are presented in Sec. IIIII The results and com- 
parisons for the forward TCs are presented in Sec. IIVI 
Section IVl provides a summarizing closure. 



II. THE MODEL AND THE BINARY 
APPROXIMATION 

We consider chains of granules all made of the same 
material of density p. When neighboring granules col- 
lide, they repel each other according to the power law 
potential 



V 



'I \r 

-TkWk - yk+i\ 



(1) 



Here yk is the displacement of granule k from its posi- 
tion at the beginning of the collision, and a is a con- 
stant determined by Young's modulus and Poisson's ra- 
tio [13, [2^. The exponent n is 5/2 for spheres (Hertz 
potential), which we use in this paper [2^. We have de- 
fined 



fc+i 



1/2 



(2) 



where is the principal radius of curvature of the sur- 
face of granule k at the point of contact. When the gran- 
ules do not overlap, the potential is zero. The equation 
of motion for the fcth granule is 



Ml 



= ar^_;^(yfc_i - y^)" 0{yk-i - yk) 

- arl{yk-yk+ir-'0{yk^yk+i), (3) 



where = (4/3)7r/9(i?J^)'^. The Hcavisidc function 9{y) 
ensures that the elastic interaction between grains only 
exists if they are in contact. Initially the granules arc 
placed along a line so that they just touch their neigh- 
bors in their equilibrium positions (no precompression), 
and all but the leftmost particle are at rest. The ini- 
tial velocity of the leftmost particle (fc = 1) is Vi (the 
impulse). Wc define the dimensionlcss quantity 



a{R[) 



I Nri+l/2 



(4) 



and the rescaled quantities Xfc, t, irik, and Rk via the 
relations 



^'1 „,l/ri 



—a't 



R'k = R[Rk, Mk = Mmik. 

Equation ([3]) can then be rewritten as 

rk-i{xk-i - XkT^^t 



(5) 



rukXk 



[Xk-l - Xk) 



rk{xk - Xfc+i)" 9{xk - Xfc+i), (6) 



where a dot denotes a derivative with respect to t, and 



rk 



2RkRk+i 
Rk + -Rfe+i 



1/2 



(7) 



The rescaled initial velocity is unity, i.e., i'i(t = 0) = 1. 
The velocity of the k-th granule in unsealed variables is 
simply Vi times its velocity in the scaled variables. 

The equations of motion can be integrated numerically, 
and we do so for our TCs. In each chain we observe that 
the initial impulse quickly settles into a propagating pulse 
whose front typically involves only three granules. In the 
backward TCs the narrow pulse is "clean" in the sense 
that granules behind it are reflected back and thus do 
not participate in forward energy propagation. In the 
forward TCs the narrow front carries a tail of moving 
granules behind it, but this tail does not affect the prop- 
agation dynamics of the narrow front. From the numer- 
ical integration results we are able to extract quantities 
such as the pulse peak location, amplitude, and velocity 
as functions of granule index and of time. These are the 
quantities to be compared with those obtained from the 
binary collision theory. 

The binary collision approximation is based on the 
assumption that the transfer of energy along the chain 
occurs via a succession of two-particle collisions. First, 
particle fc = 1 of unit velocity collides with initially sta- 
tionary particle fc = 2, which then acquires a velocity V2 
and collides with stationary particle A; = 3, and so on. 
The velocities after each collision can easily be obtained 
from conservation of energy and momentum. The veloc- 
ity Vk+i of granule fc -I- 1 after the collision with granule fc 
is always positive, so this particle goes on to collide with 
the next one. On the other hand, the velocity of granule 
fc after the collision may be positive or negative depend- 
ing on the direction of tapering of the chain (it would be 
zero for a monodispcrse chain). From these velocities, a 
number of other results follow which we outline here, but 
implement for particular chains in later sections. 

After the collision of granules fc and fc -I- 1, the latter 
emerges with velocity 



Vk+l 



2vk 



nik 



(8) 



This result can be implemented recursively to obtain 



k-l 



Vk 



Hi- 



ruk' 



(9) 



k' = l 



rrik' 



To actually evaluate this product wc need to implement 
a tapering protocol, cf. below. 

We wish to use these results to calculate the time it 
takes the pulse to move along the chain and, more specif- 
ically, the time it takes the pulse to arrive at the fcth 
granule. For this purpose, it is convenient to introduce 
the difference variable 



Zk ^ Xk - Xk+l- 



(10) 
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The equation of motion for the difference variable is ob- 
tained by subtracting the equations of motion of the two 
granules during a collision, cf. Eq. ([6]), 



Xk+l 



TTLk 



{xk - Xk+l) 



ri-1 



mk+i 



{xk - Xk+lY' 



which directly leads to 

^J.kZk 



-rk4~\ 



(11) 



(12) 



Here /i^ = mkmk+i/{mk + is the reduced mass. 

Equation p2|) is the equation of motion of a particle of 
mass /ife in the potential (rk/n)zj^ defined for Zk > 0. 
The initial condition Zk (0) is simply Vk since the velocity 
of granule fc + 1 is zero before the collision. Therefore, 
from Eq. ([9]) we can write 



fc-i 



ifc(o) = n 



=1 1 



mk'+i 

ITlk' 



The energy conservation condition 



leads to 



nfik 



1/2 



(13) 



(14) 



(15) 



We say that the pulse arrives at granule k when the 
velocity of granule k surpasses that of granule k—1, and 
that it moves on to granule fc + 1 when the velocity of the 
(fc + l)st granule surpasses that of the fcth granule. The 
residence time Tk on granule k is the time that granule 
fc takes to transfer the pulse from fc — 1 to fc + 1, and is 
given by 



Tk = 



dzk 

Zk 



dzk 



1/2 ' 



(16) 



where is the compression when the velocities of par- 
ticles fc and fc + 1 are equal [which is also the maximum 
compression, cf. Eq. (fT4)) ]. 
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The integral can be performed exactly to yield 

l/n 

/ / /, / 1. 1^ \ 

k 



- f nf^k\"''"' ,»M(-l+2/n) r(l + l/7l) 

H^j [^^(«)]' r(l/n+l/2)- 



Finally, the total time taken by the pulse to pass the fcth 
granule is obtained by summing Tk, 



(19) 



For a monodisperse chain ^k = 1/2, rk = 1 and ifcfO) = 
1, and Eq. pB]) reduces to the result obtained in Q. To 
evaluate these times explicitly for tapered chains we need 
to specify a tapering protocol. 



III. BACKWARD TAPERED CHAINS 

We proceed to implement these results for various ta- 
pered chains of spherical granules {n = 5/2), starting 
with backward tapering in this section. Two forms of 
tapering are considered. 



A. Linearly tapered chains 

In backward linearly TCs the radii of the granules grow 
linearly as Rk = 1 + S{k — 1), where 5 > is a fixed 
parameter. The ratios of the radii and of the masses of 
two successive granules then are 



Rk 



Rk+1 

rrik 
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Sk 
S 



ruk+i \ 1 + Sk 
The ratio of the masses can be substituted into Eq. 

fe-i 



Vk = 



n 



k'=i 1 



1 - 



s 



1 + Sk' 
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(20) 



(21) 



A large-fc limit can be evaluated by calculating the 
derivative of Invk with respect to fc, retaining leading 
terms in a Taylor expansion when both fc and Sk ^ 1, 
and exponentiating. This leads to the asymptotic scaling 
relation 



Vk 



fc-3/2. 



(22) 



The amplitude of the velocity pulse decreases with in- 
creasing grain number, as would be expected since the 
granules grow in size and mass. 

The associated parameters for a backward linearly TC 
are given by 



fik = {l + Skf 



1 
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1 + Sk 



Tk = {l + Sk) 
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V 1 + Sk 



(23) 



(24) 



Using Eqs. ^ and ^ in Eq. we find that 

in the large fc limit Tk varies with k according to 



Tk ~ fci3/io. 



(25) 
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The total time t taken by the pulse to reach the kth 
granule is obtained from Eq. p^ . For large k, t varies 



(26) 



This growth with k again reflects the slowing down of 
the pulse as the granules become more massive. Using 
Eq. ((26l) in Eq. (|22p . the decay in time of the pulse am- 
plitude is found to be 



(27) 



Finally, in order to compute the speed of the pulse in 
space, we note that at < = the position of each granule 
in the uncompressed chain is 



Xo{k) = Sk^ + 2k{l - 5) + 5 - 2, 



(28) 



and when the pulse has moved to the fc-th granule its 
position is given by 



X{k,t)^Xo{k)+Xk{t). 



(29) 



For large k values, the static contribution Xo(fc) domi- 
nates and the position of the pulse is 



X{t) ^20/23^ 



(30) 



The speed c{t) of the pulse therefore varies in time as 

^-3/20, (31) 



/ N dX 



Thus, the binary collision approximation in the back- 
ward linearly TC leads to a power law dependence in 
both fc, Eq. (1261), and t, Eq. ((271), for the pulse ampli- 
tude, for the position ([50]). and for the speed (PT|) . Note 
that these asymptotic decays are independent of the ta- 
pering parameter S. 

The comparison of the predictions of the binary col- 
lision approximation with the numerical integration of 
Eq. ^ is carried out as follows. First, consider the de- 
cay of the amplitude of the velocity pulse as a function of 
k. In Fig.[l]we show Vk vs k from numerical integration of 
the equations of motion (solid circles) and as calculated 
from the binary collision theory, Eq. (PT|) (open circles). 
The comparison is made for various values of 5, as indi- 
cated in the caption. Note that in some regimes of this 
figure the decay of the pulse amplitude is S'-dependent, 
as found in Eq. ()2ip and also from the numerical inte- 
gration; the 5-indcpendcnt power law decay (|22p is only 
valid when Sk ':S> 1- Figure [T] reveals the single weak- 
est aspect of the theory in that it shows substantial dif- 
ferences in the absolute magnitude of the pulse ampli- 
tude between the theoretical and numerical integration 
results, but the rate of decay of Vk is captured very well, 
within w 1% — 2%, as emphasized by the inset. The inset 
also shows the approach to the ^'-independent power 3/2 
for large fc, as predicted in Eq. 




100 



FIG. 1: The numerical integration data (solid circles) is com- 
pared with the binary collision approximation (open circles) 
for the decay of the velocity pulse amplitude with k for back- 
ward hnearly TCs with with S = 0.2, 0.4, 0.6, and 0.9, from 
top to bottom. The inset shows the decay exponent, whose 
theoretical value approaches 3/2. The numerical integration 
results also approach an S-independent value that is slightly 
larger than but close to the theoretical value. 



The decay of the velocity pulse in time as obtained 
from numerical integration is shown in Fig.O The power 
law decay (|27p predicted from the binary theory is in ex- 
cellent agreement with the numerical data. In the inset 
we show the power law exponent for different S values. 
The long-time value of the exponent 15/23 = 0.652 . . . 
obtained from the theory is in very good agreement 
(within 2% — 3%) with the observed data. Furthermore, 
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FIG. 2: The decay of the velocity pulse amplitude with time 
for backward linearly TCs with S — 0.4 to 0.9 in steps of 0.1 
from top to bottom, as obtained from numerical integration. 
The inset shows the rate of decay for different S values. The 
long-time binary prediction is 15/23 = 0.652 . . ., within 2 — 3% 
of the numerical integration results. 

gratifying agreement is found if we compare the time Tk 
obtained from Eq. (|18|) with numerical integration re- 
sults, as seen in Fig. [31 The time t that the pulse takes 
to pass by granule k is also very well reproduced by the 



5 





Ann 




500 




400 


Sh 


300 




200 




100 









100 



FIG. 3: Residence time of the pulse on granule k for backward 
linearly TCs. Comparison between the binary collision ap- 
proximation (open circles, barely visible because of the agree- 
ment between theory and numerical integration) and direct 
numerical integration of the equations of motion (solid cir- 
cles). From bottom to top, S — 0.1 to 0.9 in steps of 0.1. 



theory, cf. Fig.[4l 
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FIG. 4; Position of the pulse in units of grains versus time 
for backward linearly TCs. Comparison between the binary 
collision approximation (open circles, barely visible because 
of the agreement between theory and numerical integration) 
and direct numerical integration of the equations of motion 
(solid circles). From top to bottom, S = 0.3 to 0.9 in steps of 
0.1. 

In Table [J the agreement between the theory and the 
integration results is made exphcit. We hst the results of 
fitting numerical integration data for much longer chains 
{N > 400) ith a power law form k = a& . The numerical 
uncertainties (standard deviations) are in the third deci- 
mal place. For small values of S the data does not reach 
the values of k that are needed for the asymptotic behav- 
ior of Eq. (HSl) to be valid {Sk > 1). For the larger 
values of S the asymptotic behavior is approached, and 
here the agreement with the predicted asymptotic value 
h = 10/23 = 0.435 . . . becomes evident. Since Eq. ^ 
is valid only in the large k regime, in these fits we have 
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theory 


numerics 


0.3 


0.439 


0.438 


0.4 


0.438 


0.437 


0.5 


0.437 


0.436 


0.6 


0.436 


0.436 


0.7 


0.436 


0.436 


0.8 


0.436 


0.435 


0.9 


0.436 


0.435 



TABLE L Backward linearly TCs: Values of b obtained from 
the binary collision approximation and from the numerical 
integration of the equations of motion in the asymptotic ex- 
pression k ~ for the position of the pulse (in units of grains) 



as a function of time, 
is 10/23 = 0.435.... 



The theoretical asymptotic value of b 



only used the information for granules 300 — 400. 



B. Exponentially tapered chains 

A backward exponentially TC is one in which the ratio 
of radii of successive granules is constant, 



Rk 



R 



k+l 



= 1-9, 



(32) 



where < g < 1. The radii of granules thus increase in 
geometric progression. The radius of the fcth granule is 
= (1 - qY~^, and the mass ratio of two consecutive 
granules is 



mk 
ruk+i 



= (1 - 1? 



(33) 



This is the case that has been considered in many previ- 
ous studies [ll,[li,[i3,[ll,[ll. This mass ratio in Eq. ^ 
leads to the velocity amplitude 



where 



1 r 
2 



l + {l-q) 
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and to the associated parameters 
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r, = (l-9/2)-i/2(l 
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(34) 



(35) 



(36) 



(37) 



Substituting Eqs. ([Ml), dM]), and (071) in Eq. (HH]) and 
following the steps implemented for the linearly TC we 
find that the time t for the pulse to reach the kth granule 
now grows exponentially with k for large k, 



t ■ 



where 



77(g) = In 



V 1-9 



(38) 



(39) 
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Combining Eq. (jM]) with Eq. ([38]) leads to the decay 
of the pulse amplitude as a function of time, 



with 



f{q) 



1 



InAiq). 



The position of the fcth granule at time t is 

- l)+Xk{t). 



(40) 



(41) 



(42) 



For large fc, we can again ignore the contribution due to 
the displacement compared to the dominant static part. 
This allows us to approximate the pulse position as 



Xik) 



2-q 



(43) 



9(1-9) 
which implies that 

X{t) - (44) 

and hence the pulse velocity in time decays as 

c(0 - (45) 

Thus, for backward exponentially TCs the pulse ampli- 
tude and position vary exponentially in k, but in time 
the variations are power laws. 




FIG. 5: The numerical data (solid circles) is compared with 
the binary collision approximation (open circles) for the de- 
cay of the velocity pulse amplitude in backward exponentially 
TCs for q = 0.01, 0.02, 0.05, and 0.08 from top to bottom. 
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0.01519 


0.01522 
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0.03076 


0.03088 
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0.04673 


0.04697 


0.04 


0.06311 


0.06350 


0.05 


0.07990 


0.08047 


0.06 


0.09711 


0.09788 


0.07 


0.11477 


0.11568 


0.08 


0.13287 


0.13378 


0.09 


0.15144 


0.15275 



TABLE II: Exponential rates of decay for the velocity pulse 
in fc-space for backward exponentially TCs as obtained from 



HI and 



the binary collision approximation [In^(g) in Eq. 
the numerical integration of Eq. (|6]). 



We compare the decay rates for each q value obtained by 
fitting the numerical data in Fig. [S] with those calculated 
from the theory. The results are shown in Table II. 

The power law decay of the pulse amplitude in time 
is shown in Fig. [Sj The rates of decay are close to those 
predicted by the theory, Eq. In Table [1111 we com- 

pare the two results for different q values. Note that the 
asymptotic expressions are valid only in the long time 
limit, a requirement that involves longer times as the 




FIG. 6: Decay of the pulse amplitude with time for backward 
exponentially TCs obtained from numerical integration for 
q = 0.01, 0.02, 0.05, and 0.08 from top to bottom. 

value of q decreases. The small differences in the results 
arise mainly because Eq. (PO]) is a long-time limit (and 
not because of the binary approximation per se) . 

Finally, the theoretical residence times Tfc and the time 
t for the pulse to pass through the fcth granule are in 
excellent agreement with the numerical integration data, 
as can be seen in Figs. [7] and[S] 



Figure [S] shows the exponential decay of the pulse am- 
plitude with fc as predicted from the theory, Eq. ([M)) . 
and as obtained from direct numerical integration of the 
equations of motion. We find a small difference between 
the amplitudes predicted by theory and those obtained 
from the numerical integration, but the predicted decay 
rate is again in excellent agreement with the actual data. 



IV. FORWARD TAPERED CHAINS 

Next we consider chains where the granules are succes- 
sively smaller starting from the granule that is the sub- 
ject of the initial impulse, that is, chains with forward 
tapering. We present our binary collision approximation 
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theory 


numerics 


0.01 


1.1605 


1.1279 


0.02 


1.1673 


1.1546 


0.03 


1.1740 


1.1744 


0.04 


1.1808 


1.1820 


0.05 


1.1877 


1.1803 


0.06 


1.1945 


1.1952 


0.07 


1.2015 


1.2049 


0.08 


1.2084 


1.2145 


0.09 


1.2154 


1.2256 



TABLE III: Asymptotic rates of decay for the pulse amplitude 
in time for backward exponentially TCs obtained from the 
binary collision approximation [f{q) in Eq. (|4ip ] and those 
obtained by fitting the straight line portions of the data in 
Fig.! 




FIG. 7: The pulse residence time Tk as a function of k for 
backward exponentially TCs is shown for q — 0.01 to 0.06 in 
steps of 0.01 (top to bottom). Theory and numerical results 
are denoted by open circles and filled circles, respectively. 



results and compare them with the results obtained from 
the numerical integration of the full equations of motion 
©• We again consider two different tapering protocols, 
the linear and the exponential. 



A. Linearly tapered chains 

In forward linearly TCs the radii of the granules de- 
crease linearly as = 1 — S{k — 1). Note that since 
the radius of the first granule is unity, this places an 
5'-dependent restriction on the chain length N (or an A^- 
dependent restriction on the tapering parameter S), 



(46) 



In our calculations we set N = 100, and accordingly re- 
strict the value of S' to be < 0.01. The ratios of the radii 
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FIG. 8: The position of the pulse in units of grains versus time 
t for backward exponentially TCs is shown for q — 0.01 to 0.06 
in steps of 0.01 (top to bottom). Theory and numerical results 
are denoted by open circles and filled circles, respectively. 



and of the masses of the fcth to (fc -I- l)st granules are 
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ruk+i V l-Sk^ 
Substituting Eq. (gll) in Eq. (P, we obtain 
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(47) 
(48) 

(49) 



When k 1 + 1/S'we follow the steps described in going 
from Eq. (pij) to Eq. to find the scaling behavior 



k 



1 + 1/5 



-3/2 



(50) 



Since k < 1+1/ S, this function grows as k increases. This 
is expected since grain mass decreases with increasing k. 
The associated parameters /Xfc and for forward linear 
TCs are obtained by replacing S with —S in Eqs. ((23l) 
and ([M]) . Using Eq. ([501) ™ Eq. (US]) and summing over 
k we obtain the scaling for the time it takes the pulse to 
pass through the fcth granule, 



t ~ C(S') 



1 



(' + *>""°('-TTI7s) 



23/10' 



where C{S) ^ Co(S') [1 + 0(5')] and 

r(7/5) 



CoiS) 



lOV^ /5 
23 U 



r(9/io) 



1 



(51) 



(52) 



Figure [9] shows the variation in the pulse amplitude 
with k as obtained from the theory (|49|) and the numer- 
ical integration. As predicted from Eq. (j50p . we observe 
a linear growth in the pulse amplitude for small k. Also, 
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20 40 60 80 100 



k 

FIG. 9: The numerical integration data (filled circles) is com- 
pared with the binary collision approximation results (open 
circles) for the velocity pulse amplitude in forward linearly 
TCs for S from 0.001 to 0.008 in steps of 0.001 (bottom to 
top). The inset shows the variation of the parameter b with 
S in the fit of Vk described in the text. 



reflection of the fact that the granular mass and size be- 
come very small and as a result the pulse passes through 
these grains very quickly. 

The change in the pulse amplitude with time is shown 
in Fig. [TlJ The theoretical result for the time-dependent 
amplitude can be obtained by combining Eq. ([5n|) with 




20 40 60 80 100 

t 



a.s k —>■ N the numerical results show that the velocity 
behaves as Vk ~ a{l — bk)~^^^ and thus diverges around 
k ^ 1/6. Interestingly, we find that 1/6 is very close to 
(1 -I- 1/5) and thus Eq. ([5D|) . which is in principle re- 
stricted to the regime fc <C 1 + 1/5, seems to work well 
all the way up to the divergence. In the inset in Fig. [5] 
we compare the values of 6 obtained from the theory, 
6 = (1 -I- 1/5)~^, and the fit to numerical data. 

In Fig. [ini we compare the results for the time t taken 
by the pulse to pass through the fcth granule. The 




t 



FIG. 11: Change in the pulse amplitude with time for forward 
linearly TCs. as obtained from numerical integration of the 
equations of motion. Results are for S from 0.001 to 0.01 in 
steps of 0.001 from bottom to top. 

Eq. dHI). We find that 

/ X -15/23 

"<"^('-C(5)) ■ <^=" 

For small t the pulse amplitude thus grows linearly with 
t, and as t — > C(s) it diverges with exponent 15/23. Since 
the absolute value of the pulse amplitude is not captured 
correctly by the binary collision approximation, the di- 
vergence time C(5) of the pulse amplitude is also differ- 
ent from the actual divergence time seen in the numerics. 
However, we find that the initially linear behavior and 
the divergence exponent 15/23 are accurately predicted 
by the theory. In Fig. [T^] we show the pulse residence 
time obtained for each granule from numerical integra- 
tion along with the corresponding values obtained from 
the binary collision approximation. 



FIG. 10: Comparison of the theory (open circles) and numer- 
ical integration results (filled circles) for the time t taken by 
the pulse in forward linearly TCs to reach the fcth granule for 
S from 0.001 to 0.01 in steps of 0.001 from right to left. 

change in t with k is described extremely accurately by 
the binary collision approximation result (|5ip valid for 
k <C 1/5. The result in Eq. ([5T|l is thus expected to be 
valid for all k provided that 5 is small (1/5 3> 1). Inter- 
estingly, the result in Eq. (|5T|) in fact works well for all 
k < N (with N = 100) in the range 0.001 < 5 < 0.01. 
The time t increases linearly with k for small k but be- 
comes almost /c-independent as k increases. This is a 



B. Exponentially tapered chains 



In this case the radii of the granules decrease as 



Rk = 



(1 



(54) 



where q > 0. The velocity of the fcth granule can be 
obtained from Eq. by replacing q by —q. Figure [T^ 
shows the comparison of the binary collision approxima- 
tion results with the numerical data. We again observe 
that although the absolute values differ by a constant 
amount, the rate of growth is correctly predicted by the 
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FIG. 12: Residence time of the pulse on each granule for for- 
ward linearly TCs for S = 0.001, 0.002, 0.003, 0.005, and 
0.009 from top to bottom. Numerical integration results: 
filled circles; binary collision approximation results: open cir- 
cles. 




q 


theory 


rmmerics 


0.01 


0.01481 


0.01471 


0.02 


0.02926 


0.02879 


0.03 


0.04336 


0.04267 


0.04 


0.05710 


0.05565 


0.05 


0.06848 


0.07051 


0.06 


0.08359 


0.07833 



TABLE IV: Comparison of the rates for the exponential 
growth of the pulse amplitude in forward exponentially TCs 
in granule space obtained from the binary collision approxi- 
mation with those obtained by fitting the data in Fig. 1131 



100 




FIG. 14: Total time taken by the pulse to pass through the 
kth granule in forward exponentially TCs for q from 0.01 to 
0.06 in steps of 0.01 (right to left). Numerical integration 
results: filled circles. Binary collision approximation results: 
open circles. 



FIG. 13: Change in the pulse amplitude as a function of k for 
backward exponentially TCs for q from 0.01 to 0.06 in steps of 
0.01 from bottom to top. Numerical integration results: filled 
circles. Binary collision approximation results: open circles. 



theory. In Table HVl wc compare the rates obtained from 
the fit to the data points in Fig. [T3] with the results of 
the theory. 

The total time taken by the pulse to reach the kth 
granule is obtained by substituting Eqs. ([M]) and (^5]) 
in Eq. (jlSp . replacing q by —q, and summing over k, cf. 
Eq. ([19]). For long times, we obtain 



t ~ 1 



,(fe-l)r,(-g) 



(55) 



where r](—q) is obtained by replacing q by —q in Eq. (|39p . 
Since ri{—q) < 0, t approaches a constant value at large 
k. In Fig.[T4]we show both the theory and the numerical 
results for the time taken by the pulse to pass through 
the fcth granule. 

The time dependence of the pulse amplitude in the 



binary collision approximation is obtained as 



vit) ~ 1 + 



M{q) 



(56) 



where 



M{q) = 



2V^ r(7/5) 5 ^n^q/2 



2/5 



5r/(-q) r(9/10) UH-(l + g)3 



, (57) 



and f{—q) is obtained by replacing q with —q in Eq. (j4ip . 
At short times v{t) increases linearly with t. Also, since 
ri{-~q) < 0, M{q) < and, as t — ^ M{q), the velocity 
diverges with the exponent f{—q). These behaviors are 
shown in Fig. [121 Again, since the absolute value of the 
pulse amplitude is not captured correctly within the bi- 
nary collision approximation, the divergence time of the 
pulse amplitude M{q) is also different from the actual 
divergence time seen in the numerics. However, the ex- 
ponent f{—q) is accurately predicted by the theory. In 
Table |V| we compare the exponent f{—q) obtained by 
fitting the numerical data in Fig. [TSj with the results of 
the theory. 



10 



10^ 




10 20 30 40 50 60 70 
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FIG. 15: Numerical integration results for the pulse amplitude 
as a function of time in forward exponentially TCs for q from 
0.01 to 0.06 in steps of 0.01 (bottom to top). 



q 


theory 


numerics 


0.01 


1.1472 


1.1361 


0.02 


1.1406 


1.1308 


0.03 


1.1341 


1.1134 


0.04 


1.1276 


1.1237 


0.05 


1.1211 


1.1022 


0.06 


1.1147 


1.1101 



TABLE V: Comparison of the power law exponent f{—q), 
Eq. (|56p . obtained from the binary collision approximation 
with the exponent extracted from numerical integration re- 
sults for forward exponentially TCs. 

V. CONCLUSIONS 

We have introduced an analytic approach to calculate 
the behavior of a propagating pulse along a variety of 
tapered chains of spherical granules. We implement a bi- 
nary collision approximation Q that supposes that the 
collision events by which a pulse propagates involve only 
two granules at a time. While the approximation overes- 
timates the pulse amplitude, it captures all other prop- 
erties remarkably well. These include the rate of decay 
of the pulse amplitude with time, the residence time of 
the pulse on each granule, and the pulse speed in units 
of granule and in space. 

We have implemented our approximation in four dif- 
ferent tapered chains: backward linearly and exponen- 
tially tapered chains (the granules increase in size), and 
forward linearly and exponentially tapered chains (the 
granules decrease in size). In backward linearly tapered 
chains, where granular size increases linearly along the 



direction of pulse propagation, an initial impulse on the 
first granule settles into a pulse whose amplitude, speed 
and width change slowly as the pulse propagates. The 
pulse speed and the pulse amplitude decay in granular 
units and also in real time, as is reasonable since the 
granules increase in size and mass along the direction 
of propagation. Although the pulse remains narrow in 
granule number, the pulse width in real space of course 
increases as the granules get larger. All these decreases 
and increases vary as power laws both in granular units 
k and in time t for large k and t, and the decay and 
growth exponents are insensitive to the value of the ta- 
pering S over the range that we have tested, 0.1 < 5 < 1. 
On the other hand, for backward exponentially tapered 
chains, where granule size increases geometrically, the 
pulse properties change exponentially in k while they ex- 
hibit power law behavior in time. In this case the rates 
arc strongly influenced by the tapering parameter q. 

In forward linearly tapered chains, where granular size 
decreases linearly, the pulse speeds up both in granule 
number units and in time, and in fact diverges. The di- 
vergence exponents are again insensitive to the tapering 
parameter S and are quantitatively reproduced by the 
theory. Similar divergences are observed in the forward 
geometrically tapered chain, but now the behavior is ex- 
ponentially sensitive to the tapering parameter q. 

In the binary collision approximation the initial veloc- 
ity for the collision between granules k and A; -I- 1 is taken 
to be the velocity of granule k at the end of the colli- 
sion between granules k — 1 and k. This estimate in turn 
leads to a pulse velocity amplitude that is higher than 
that obtained from the numerical integration of the full 
equations of motion, but all other pulse characteristics 
are quantitatively reproduced in all our tapered chains 
by this analytic approximation. Our next challenge is to 
generalize this approach to chains with other mass vari- 
ation profiles and even to random chains. This work is 
in progress. 
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